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We present a new indirect method to measure the quantum state of a single mode of the elec- 
tromagnètic field in a cavity. Our proposal combines the idea of (endoscopic) probing and that of 
tomography in the sense that the signal field is coupled via a quantum-non-demolition Hamiltonian 
to a meter field on which then quantum state tomography is performed using balanced homodyne 
detection. This technique provides full information about the signal state. We also discuss the 
influence of the measurement of the meter on the signal field. 

PACS numbers: 03.65.Bz, 42.50.Dv 



I. INTRODUCTION: STATE MEASUREMENT 

The question of how to measure the quantum state 
of a single mode of the electromagnètic field in a cavity 
has recently attracted a great deal of attention M, as 
it determines the feasibility of the measurement of the 
state of a genuine quantum system Several propos- 
als have been made in the last few years to answer this 
question. Among them, quantum-state endoscopy f| and 
quantum-optical homodyne state-tomography [0] are two 
notable examples. The former proposal makes use of 
a beam of two-level àtoms — sent with controlled speed 
through the cavity — to infer the properties of the field 
state inside the cavity. On the other hand, optical homo- 
dyne tomography is a method for obtaining the Wigncr 
function (or, more generally, the matrix elements of the 
density operator in some representation) of the electro- 
magnètic field. It therefore consists of an ensemble of re- 
peated measurements of one quadrature operator for dif- 
ferent phases relative to the local oscillator of the homo- 
dyne detector. The major drawback of the first method 
is the low detection efhciency for àtoms, whereas in the 
second one one has to couple the field out of the res- 
onator. 

In the present paper we propose to couple the field 
via a quantum-non-demolition (QND) interaction || to 
a meter field on which we then perform tomography using 
a balanced homodyne detector. In this way we combine 
the idea of probing, that is doing endoscopy on the field 
without taking it out of the cavity, and the tool of tomog- 
raphy and arrive at the method of endoscopic quantum 
state tomography. In contrast to the method of quantum 
state tomography pj based on homodyne detection, the 
present technique does not couple the signal field out of 
the resonator. 



II. OUR MODEL 

Let us assume that the electromagnètic field we want 
to probe (the signal mode) is in a pure quantum state and 
neglect dissipation. We note, however, that our method 
applies as well to a field described by a density operator. 
In order to measure the signal field we couple it to a meter 
field. Both fields are then coupled to a pump field. This 
coupling leads us to a quantum-non-demolition Hamil- 
tonian describing the interaction between the signal and 
the meter mode. 

Our model starts from the Hamiltonian 



H = ihx[a. 
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where à s (à\), àmia^), and à p (àj ) ) represent the annihila- 
tion (creation) operators of the signal, meter, and pump 
field, respectively. The parameters x an d cr measure the 
coupling between the three fields, and the meter and sig- 
nal field, respectively. When the pump field is highly 
excited we can describe it by a coherent state of ampli- 
tude a and phase 2<p, that is 



(2) 



If we substitute the coherent state approximation, 
Eq. (||), into the Hamiltonian, Eq. (Q), after some àlgebra 
we obtain 



H = 2haX s ((f> + tt/2) ■ X m (<j>) , 



(3) 



where we have arranged the strength a of the pump field 
such that xoí — <?, and the quadrature operators are given 
by 



y/2 S 



(4) 



of the signal (j = s) and the meter (J = m) mode at the 
angle 9. 
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It is particularly interesting to note that due to the 
special choice x a — a we have been able to obtain an in- 
teraction between the signal and the meter which couples 
the quadrature operator X m (4>) of the meter at phase an- 
gle 4> to the out-of-phase quadrature operator X s (0+7r/2) 
of the signal . The present Hamiltonian is a genuine 
QND Hamiltonian. In the next sections we shall see how 
such a Hamiltonian can be used to measure the quantum 
state of the signal field. We note that according to the 
QND Hamiltonian Eq. (0) a (homodyne) measurement 
of the meter at a fixed phase <j) of the pump field yields 
information about the signal in the orthogonal quadra- 
ture. By varying the phase <p of the pump field we can 
probe all quadratures of the signal. 



III. ENTANGLED STATE 

The aim of the present section is to calculate the 
combined state l^) of signal and meter achieved after 
some interaction time according to the QND interaction 
Hamiltonian, Eq. (Q). 

When we couple for an interaction time t the signal 
and meter mode prepared initially in the states \ip s ) and 
\ipm) we find the quantum state 

|*(í)) = exp(-iÈt/h)\il> m )\il),) 

= exp[-2iatX s (<f> + 7r/2)X m (<£)]|V> m )|VO (5) 

for the combined system. 

We expand the initial signal state in quadrature states 
|X s (0 + 7r/2)>, 

OC 

|^)= [ dX s A(X s ;<t> + n/2)\X s (cf> + Tr/2)) . (6) 

— oo 

We stress that this representation and, in particular, the 
wave function tjj s (X s ; (f>+ir/2) = (X s ((f)+iT/2)\tp s ) depend 
crucially on the angle 9 S . 

Combining Eqs. (^|) and (Q), we may rewrite the com- 
bined state as 

oo 

|*(t)> = / dX s iP s (X s ;4> + Tr/2)\X s (4> + Tr/2)) 



x exp[-2iatX s X m ((f))]\>(p n 



(7) 



Expanding \ip m ) in quadrature states \X m (9)} of the 
meter at the angle 9, that is 



\Í> m )= / dX m ^ m {X m -9)\X m (9)) 



(8) 



where , m (X m ;0) = (X m (9)\ip m } denotes the wave func- 
tion of the meter state at the angle 9, it is straightforward 
to find S 



exp [-i(2atX a )X m (<f>)]\ij) m ) 

dX rn exp[-ij(X s , X m ; 9 - </>)] (9) 
x i> m [X m - 2atX s sm(9 - 0); 9}\X rn (9)) , 

where 

j(X s ,X m ; 9-4>) = {atX s f sin[2(0 - 4>)] 

+2<rtX s X m cos(9 - (/>) . (10) 

Hence the combined quantum state reads 

oc oo 

|*(í))= j dX s J dX m i> s (X s ;<p + n/2) 

— oo — OO 

xip m [X m - 2<jtX s sm(9 -4>); 6] 

x exp[-Í7(X s , X m ; 9 - cj>)] 

x\X s (<t> + n/2))\X m (9)) . (11) 

We note that due to the coupling between the meter and 
the signal via the Hamiltonian Eq. (^|), the meter wave 
function %j) m (X m \9) at the angle 9 gets shifted by an 
amount SX m = 2atX^ sm{9 — ó). 



IV. EFFECT OF THE METER MEASUREMENT 
ON THE SIGNAL STATE 

In the present section we shall show how a measure- 
ment of the meter influences the state of the signal. Let 
us first consider an arbitrary quadrature state of phase 
angle 9. 



According to Eq. (11) the conditioned state 
1 



y/W(X m ) 



(X m (8)\*(t)) 



(12) 



of the signal given that our quadrature measurement at 
angle 9 has provided the value X m reads 



|^ c) )= / dX s 4, s {X s ;^ + -K/2)f{X s \X m )\X s {cjy + TT/2)) , 

— OO 

(13) 

where the "filter function" / is given by 
1 



f{Xg\X m ) — 



ip m [X m -2atX s sm(9-^y,e] 

(14) 



VW{X. 
x exp[-i^(X s , X m ; 6 - cj>)) 



The normalization condition dircctly yields he probabil- 
ity W(X m ) of finding the meter variable X m , that is 
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oc oc 



W{X V 



dX s |^ s (X s ;0 + 7r/2)| 2 



x \^ m [X m ~2atX s ún{9-<f ) )-9]\ 2 . (15) 



Equation (13) clcarly shows how the measurement of thc 



mcter influences the quantum state of the signal: The 
filter function determined by the wave function of the 
meter selects those parts of the signal wave function that 
are entangled with the corresponding parts in the meter. 
To study this in more detail we now calculate the Wigner 
function 

oo 

W^(X s ,P s \X m ) = -L / dYe lP ° Y (X s -Y/2\^) 
2ir J 

— oo 

x(4 c ï\X s + Y/2) (16) 

of the signal state conditioned on the measured meter 
value X. m . Substituting the state l^s ), Eq. dÏ3|), into 
this expression we arrive at 

oo 

WÍ c \X s ,P s \X m ) = ±- J dYe iP ° Y i> s {X s -Y/2) 

— oc 

x r s {X s + Y/2)f(X s - Y/2\X m )f*(X s + Y/2\X m ) . (17) 
The integral may be expressed as the convolution 
W^(X s ,P s \X m ) = 

oo 

í dP'W s (X s ,P s -P')W f (X s ,P'\X m ) (18) 

— oo 

between the Wigner function of the original signal state 
W S {X S ,P S ) = 

OC 

i- j dYe* p ° Y i> s (X s -Y/2)r s (X s + Y/2) , (19) 

— OO 

and the Wigner function 
W f {X s ,P s \X m ) = 



2tt 



dye tP ° Y f(X s - Y/2\X m )f*(X s + Y/2\X rn ) (20) 



of the "filter" provided by the measurement on the meter. 



V. SPECIAL EXAMPLES 



A. In phase measurement 



If we take the angle 9 equal to 4>, the state | , Eq. (11) 
of the complete system reduces to 



|*(í))= / dX s J dX m i/j s (X s ;4> + Tr/2)i> m (X m ;<t>) 

— OO —OC' 

xexp(-i2atX a X m )\X s (4> + ir/2))\X m (4>)) . (21) 

In this case the meter wave function is not shifted. Never- 
theless, the states of signal and meter are still entangled. 
Since the shift SX m vanishes, the probability 

oo 

W(X m ) = \^ m \ 2 í dX s \MX s )\ 2 = \^rn\ 2 , (22) 



of finding the meter variable X m following from Eq. ( p^| ) 
for 9 — 4> is identical to the initial probability of the 
meter, that is 



W(X m ) = \i/> m (X m ) | 



(23) 



Hence, up to an overall phase fi m determined by the 
meter wave function ip(X m ) = \ip(X m )\exp[ifj,(X m )], 
we find from Eq. ( |Ï4| ) thc filter function f(X s \X m ) = 
exp(— i2atX s X m ), and from Eq. ( |Ï3[ ) the conditioned sig- 
nal state 

oo 

|^ c) ) = J dX s ^ s (X s )exp(-i2atX s X m )\X s (^ + n/2)) . 

— OO 

(24) 

Note that the measurement of the meter has indeed 
changed the state of the system but did not alter the 
probability 

W(X S ) = |(X S |^ C ))| 2 = \^ s (X s )exp(~ ï 2 ( jtX s X m )\ 2 
= \M*s)\ 2 (25) 

of finding the signal variable X s . Hence, the measure- 
ment has left untouched the shape of the original state 
(in the Wigner function representation) but has moved 
it along the momentum axis by an amount of 2atX m . 
Consequently, the measurement did not change the prob- 
ability distribution in the conjugate variable, namely the 
X s variable. We note, however, that in this way we can- 
not gain information about the signal since according to 
Eqs. @ and (§jj) the probability distribution W(X m ) 
of measuring the variable X m is identical to the original 
distribution. 

This finding is actually a rather general result. In fact, 
it can be rigorously shown Q that a (QND) measurement 
which does not change the probability density of the ob- 
servable which is being measured on a single quantum 
system gives no information about the measured observ- 
able. 
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B. Out of phase measurement 



B. Tomographic measurements 



Turning now to the case of 8 — <fi + tt/2, we see that 
the shift 8X m = 2atX s in the meter wave function is 
maximal and according to Eq. ( |ïo| ) the phase 7 vanishes. 
Hence, the combined state 



dX m ip a (X s ; + tt/2) 



(26) 



In the present section we show that it is possible to 
perform tomography on the meter mode to obtain infor- 
mation about the signal state. To this end, we rewrite 
Eq. @ 



W(X m )= / dX 8 \^ 8 {X s )\ 2 \^ m {X m -2<jtX s )\ 2 , (31) 



x*l>m{X m ~ 2atX s ; <j> + tt/2) \X s {cj> + tt/2)) \X m ((f> + tt/2)) 

is an entangled state in which the entanglement between 
the meter and signal is due to the s hift o f the meter. 
In contrast to the discussion of Sec. VA we can now 



deduce properties of the signal from the shift of the meter 
wave function. Unfortunately, we cannot simultaneously 
keep the probability distribution W(X S ) = \ifj s (X s )\ 2 of 
the original signal state invariant, in accordance with the 
discussion at the end of Sec. |VA|. Indeed, we find from 
Eqs. @ or 



the conditional state 



which gives the marginal distribution of the meter (prob- 
ability distribution of the results of the measurements of 
X m ) in the case of out of phase measurements. Let us as- 
sume that the meter wave function is extremely narrow, 
that is the meter is initially in a highly squeezed state, for 
example a squeezed vacuum |0,r), where r is the (real) 
squeezing parameter. Then, according to Eq. (|3Ï]), the 
marginal distribution W(X m ) is given by a convolution 
of the modulus square of the signal wave function with a 
narrow Gaussian 



(c)\ 



1 



W(X m ) 



dX s ip s (X s )ip m (X m - 2atX s ) \X S ) \^ m (X m - 2atX s )\< 



1 



y / 7rcoshr(l — tanhr) 



(32) 



(27) 

of the system given the meter measurement at phase 4> + 
7r/2 has provided the value X m . The probability 

00 

W{X m )= J dX s \MXs)\ 2 \^ m (X m -2atX s )\ 2 (28) 

—00 

of finding the meter value X m following from Eq. (|ï^) is 
now a convolution of the system and the meter function. 

VI. SPECIAL MEASUREMENTS 
A. Weak measurements 

If tpm is broad compared to ip s we can evaluate ip m 
at some characteristic value of X s , such as (X s ). As 
a consequence the conditional state, Eq. (|27j), is simply 
given by 



x exp 



1 + tanh j 
1 — tanh ? 



(X m — 2atX s y 



Now, if the squeezing parameter r is large enough, the 
Gaussian (^) approaches a delta function in the meter 
and signal variables 



\il) m {X m -2atX s )\ 2 - 

and Eq. (|3^) reduces to 

00 

1 



1 



2at\ 



X, 
~2at 



S \ X s - TTl) , (33) 



W{X V 



\2at\ 
1 

2M 



Í's 



dX s \ÍJ s (X s )\ 2 S[X, 
X 



2at 



2at 



2at \ 2at 



(34a) 
(34b) 



dX s ip s {Xs)\X s 



The probability 

W(X m ) = \^ m (X m + 2at(X s ))\ 2 



(29) 



(30) 



reduces to the original meter probability shifted by an 
amount 2at(X s ). Hence, when this shift 2at(X s ) is larger 
than the width of W m (X m ) = \4>m(X m )\ 2 , we can learn 
about (X s ). As seen from Eq. (|29|), in this case the state 
of the signal mode does not change appreciably. 



Hence, by measuring the probability distribution W(X m ) 
of the outeomes of the meter variable X m (for example 
via balanced homodyne detection performed on the me- 
ter field) we indirectly obtain the probability distribution 
W(X S ), up to a rescaling given by the factor 2at. How- 
ever, from Eq. ( |27j ) it is clear that in this case the signal 
wave function is changed, and therefore we need to pre- 
pare the signal field again in the same state after each 
measurement. This is what is usually done in quantum 
optical tomography [Q. 

The advantage of the present scheme is that we per- 
form an indirect measurement: We do not detect the 
signal mode outside the cavity (that is, we do not have 
to take the signal field outside the cavity), but we cou- 
ple it to a meter field which is successively detected, 
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thus overcoming the smearing effect introduced by the 
direct detection of the signal B. Moreover, there is no 
need of a smoothing procedure, since we are interested 
in the marginal probability distribution W(X S ) which is 
directly related to W{X m ) through Eq. (|34|). In order to 
probe the full state of the signal field, however, we would 
need to measure the probability distribution W{X m ) for 
various vàlues of the phase . 

VII. CONCLUSIONS 

To summarize, we have suggested a method to mea- 
sure the quadrature probability distribution (or, more 
generally, the full quantum state) of a single mode of the 
electromagnètic field inside a cavity. It is based on indi- 
rect homodyne measurements performed on a meter field 
which is coupled to the signal field via a QND interaction 
Hamiltonian. We have named this procedure "endoscopic 
tomography" because (i) it does not require (in contrast 
to Ref. 0) to take the field out of the cavity, just as in 
"quantum state endoscopy" ||, where a beam of two- 
level àtoms is used as a probe; (ii) tomographic mea- 
surements performed (by balanced homodyne detection) 
on the meter mode allow us to reconstruct the marginal 
probability distribution of the signal variable or even the 
full quantum state. 
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